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ABSTRACT 



Context. Episodic flaring activity is a common feature of X-ray pulsars in HMXBs. In some Be/X-ray binaries flares were observed 
in quiescence or prior to outbursts. EXO 2030+375 is a Be/X-ray binary showing "normal" outbursts almost every ~46 days, near 
periastron passage of the orbital revolution. Some of these outbursts were occasionally monitored with the INTEGRAL observatory. 
Aims. The INTEGRAL data revealed strong quasi-periodic flaring activity during the rising part of one of the system's outburst. Such 
activity has previously been observed in EXO 2030+375 only once, in 1985 with EXOSAT. (Some indications of single flares have 
also been observed with other satellites.) 

Methods. We present the analysis of the flaring behavior of the source based on INTEGRAL data and compare it with the flares 
observed in EXO 2030+375 in 1985. 

Results. Based on the observational properties of the flares, we argue that the instability at the inner edge of the accretion disk is the 
most probable cause of the flaring activity. 
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1. Introduction 

Accreting X-ray pulsars, mostly residing in high-mass X-ray bi- 
naries (HMXBs), often show abrupt increases in their X-ray lu- 
minosity lasting from a few tens of seconds to several hours 
- X-ray flares. Flaring activity is often observed on top of a 
slower flux variation related to X-ray outbursts or super-orbital 
modulation. Amo ng the best-known "flaring" X-ray pulsars 
are LMCX-4 (e.g. Moonetal. 2003b, end references therein), 
SMCX-1 dAn gelini et al. 1991; Moonetal. 2003a), VelaX-1 
jKrevkenbohm et al.1 120081: iFiirst et alj I2O1O0 . X-ray binaries 
with a donor star of Be (or Oe) spectral type (Be/X-ray binaries 
or BeXRB) are cur rently the most numerous class of HMXBs 
with X-ray pulsars ( Liu et al. 2006), although another class of 
HMXB pulsars, the supergiant fa st X-ray transients (SFXT), is 
rapidly growing (e.g. lSidohl20Tlh . BeXRBs are characterized by 
periodic or sporadic X-ray outbursts lasting from several days to 
several weeks when the neutron star accretes matter fr om the 
equatorial disk around the donor star (see e.g. Reig 2011, for a 
recent revie w). These systems are also known to show occasional 
X-ray flares dFinger et al.ll"999l:lReig et alj2008l:ICabaUero et alj 



l2008LlPostnov et al.ll2008 

Probably the most remarkable flaring activity among 
BeXRBs was exhibited by the 42 s pulsar EXO 2030+375. A se- 
ries of six roughly equidistant flares with a mean period of -four 
hours was observed with the EXOSAT satel lite a few months af - 
ter the giant outburst of the source in 1985 ('Parmaret al."l989|). 
EXO 2030+375 is one of the most regularly monitored BeXRBs. 
In addition to two "major" outbursts in 1985 and 2006 with a 
peak X-ray luminosity of Lx > 10""^ erg s"', assuming a distance 



of 7 kpc (IWilson et al.ll2002h . the source exhibited less powerful 
"normal" outbursts with Lx ~ 10^^ ergs '' almost every or bit. 
The orbital period of the system is ~46 d d Wilson et al.ll2008h . 

Some flaring of EXO 2030+375 has been reported since the 
EXOSAT observations, but it has generally a ppeared to be a sin- 
gle flare during the rise of normal outbursts dReig & CodI 19981 : 
ICamero Arranz et"al]|2005h . Flares were also apparently present 
during and between the normal outbursts sho rtly after the 2006 
major outburst (Fig. 3 in I Wilson et al.ll2008b . But the short du- 
ration and broad spacing of the observations make it unclear 
whether this is the same quasi-periodic phenomenon that was 
seen with EXOSAT. In this Letter we present the INTEGRAL 
observations of the source that, for the first time since the 
EXOSAT observations, reveal strong quasi-periodic flaring be- 
havior. 



2. Observations and data processing 

Since its major outburst in 2006, EXO 2030+375 has repeat- 
edly appeared in the field of view of In ternational Gamma- 
Ray Astrophysics Laboratory (INTEGRAL. IWinkler etani2003h 
during the observational program concentrated on the "Cygnus 
Region". The INTEGRAL observatory has three main X-ray 
instruments: (i) the imager IBIS sensitive from ~20keV to a 
few MeV with hig h spatial and moderate spectral resolution 
dUbertini et al .' 2003); (ii) the spectrometer SPI, which is sen- 
sitive in the roughly same energy range as IBIS, but with much 
higher spectral resolution an d substantially lower imaging capa- 
bility dVedrenne et alJl2003h : and (iii) the X-ray monitor JEM- 
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Fig.l. INTEGRAL observations of EXO 2030+375 (vertical 
bars) superposed on the Swift/BAT light curve of the source 
since its major outburst in June-October 2006. 



X with moderate spectr al and spatial reso lution, operating be- 
tween ~3 and ~35keV (iLund et al.ll200l . INTEGRAL obser- 
vations normally consist of a series of pointings called Science 
Windows, 2 to 4 ksec each. 

Because the times of the observations were not optimized 
for EXO 2030+375, the INTEGRAL coverage of the source's 
normal outbursts is rather sparse. The Swift/BAT light curve of 
the pulsafl with the indicated INTEGRAL observations is shown 
in Fig.[T] 

For our analysis we used the ISGRI de tector layer of IBIS 
sensitive in the 20-300 keV energy range dLebrun et al.ll2003h 
and JEM-X. Owing to limited count-rate statistics, no additional 
information could be gained from SPI data. The standard data 
processing was performed with version 9 of the Offline Science 
Analysis (OSA) sof tware provided by the IN TEGRAL Science 
Data Centre (ISDC. ICourvoisier et aLll2003l) . We performed an 
additional gain correction of the ISGRI energy scale based on 
the background Tungsten spectral lines. 

To search for the flaring activity, we examined the entire 
IBIS/ISGRI light curve by combining the publicly available 
ISGRI data products in the HEAVENS data bas^| with the re- 
sults of our own analysis. We did not find any clear evidence 
of flares in all the data except for the latest INTEGRAL obser- 
vations of the source in November and December 2010 (MJD 
~55520-55540), partially covering a normal outburst. The ris- 
ing part of the outburst is shown in Fig.|2l The upper panel of 
Fig. [3] shows the entire outburst as observed with INTEGRAL. 
In this work we concentrate on the analysis of the source's flar- 
ing behavior during this outburst. 



3. Timing analysis 

As can be seen in Fig.|2l the flux in the rising part of the out- 
burst experiences quasi-periodic oscillations/flares that cease as 
the averaged flux increases. One can identify at least five subse- 
quent flares with a mean period of ~0.3 days (~7 hours). 

To characterize the level of the flux variability during the out- 
burst we calculated the normalized excess variance in the light 
curve in relatively broad time intervals. We defined four equal 



' We used the Swift/BAT transient monitor resuUs provided by the 
Swift/BAT team 
^ http://www.isdc.unige.ch/heavens/heavens 
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Fig. 2. The INTEGRAL/ISGRI light curve of EXO 2030+375 
during the rise of its normal outburst in November-December 
2010. The flaring activity is clearly seen. The light curve of the 
entire outburst is shown in Fig. [3] One Crab corresponds to ~260 
cts/s in the specified energy range. 
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Fig. 3. Top: ISGRI fight curve of the EXO 2030+375 out- 
burst approximated with a polynomial function (solid curve) 
which represents the averaged evolution of the flux. Bottom: 
Normalized excess variance of the source flux with respect to 
the polynomial function. 



adjacent intervals in the rising part and five equal adjacent inter- 
vals in the decay. The normalized excess variance is often used 
as a simple mea sure of the intrinsic variability amplitude in fight 
curves, (see e.g. lNandra et al.iri997h : 



1 



' NXS 



Nif}'- 



r/ r /■aver\2 2n 



(1) 



Here is the number of data points in the corresponding time 
interval, fi is the flux of the individual data points, cr, - their 
uncertainty, fl^"" is the smoothed evolution of the flux obtained 
by a polynomial fit to the light curve (upper panel of Fig.O, 
and {/) is the mean value of the flux within the interval. The 
normalized variance cr^xs ^'^ '^^^ '-^^^ represents the amplitude 
of intrinsic flux variations superimposed on the smoothed flux 
development. The term cr, under the summation (E) eliminates 
the contribution of the Poisson noise. The bottom panel of Fig. [3] 
represents the normalized variance as a function of time during 
the outburst of EXO 2030+375. It can be seen that the amplitude 
of the variability is high at the rising phase (corresponding to the 
flaring episode). Then it decreases towards the maximum of the 
outburst and remains low during the decay. 
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Fig. 4. Top: The Lomb-Scargle periodogram of the "flaring" 
part of EXO 2030+375 light curve (between MJD 55518.5 and 
55521.5). The peak around ~0.3 d is clearly seen. Bottom: The 
"flaring" part of the light curve folded with best period found 
from the periodogram (0.293 d). 



To study periodicity of the flares, we performed a formal 
period search in the rising part of the outburst (between MJD 
55518.5 and 55521.5) using the Lomb-Scargle periodogram. 
The results are presented in Fig. ID The periodogram shown in 
the top panel indicates a clear peak around ~0.3 days (~7 hours). 
The averaged profile obtained by folding the light curve with this 
period is shown in the bottom panel. The profile shape is asym- 
metric and characterized by a steep rise and a slower decay. 

We used the INTEGRAL data to study the pulse period be- 
havior and pulse profiles during the outburst. The photon ar- 
rival times were converted to the reference frame of the solar 
system barycenter and corrected for the binary orbital motion 
using the ephemeris bv IWilsorTe t al. (2008). Using the pulse- 
phase-connection technique (e.g. lStaubert et al. 2009), we found 
the pulse period P = 41. 3 15 16(2) s at the epoch ro(MJD)= 
55526.056994 and the period derivative P = -1.9(1) x 10"'' s/s, 
that indicates significant spin-up. We used the measured pulse 
ephemeris to construct and study the pulse profiles of the source. 
We could not find any difference between the profiles obtained 
during the rise and decay of the outburst. Figure|5]shows the pro- 
files accumulated during the entire outburst. 

4. Spectral analysis 

For the spectral analysis we used JEM-X data between 3.5 and 
35keV, and ISGRI data between 20 and 80keV. We added sys- 
tematic uncertainties at the level of 2% to the JEM-X spectra 
and 1 % to the ISGRI spectra based on the recommendations of 
flie instrument teams and the Crab observations. The spectrum 
of the source during the outburst (accumulated between MJD 
55518.5 and 55533.0) was modeled with the cutoff'-powerlaw 
model (F{E) oc E^^ x exp [E/Efou], where E is the photon en- 
ergy, r and iifoid are the photon index and the folding energy, 
respectively) modified at lower energies by photoelectric absorp- 
tion. 

The best-fit parameters are F = 1.6(1), /Sfoid = 30(2) keV, the 
absorption column density «H - ll(l)x 10^^ hydrogen atoms per 
cm^. The uncertainties in parentheses refer to the last digit and 
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Fig. 5. Pulse profiles of EXO 2030+375 obtained with ISGRI 
(top) and JEM-X (bottom) instruments during the outburst (MJD 
55518.5-55533.0). 
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Fig. 6. Photon index F as a function of flux during the rising (flar- 
ing) part of the outburst (left) and the decay phase (right). One 
Crab corresponds to ~260 cts/s in the specified energy range. 



are at the 90% confidence level. The value of ne is substantially 
higher than measured in previous observations of the source, in- 
cluding older INTEGRAL measurements (1-3x10^' cm~^, e.g. 
iKlochkov et al.1 120071: IWilson et al] |2008|) . We note, however, 
that studying riu with JEM-X is generally problematic as the data 
below ~3 keV are not available. The significance of the mea- 
sured increase in absorption is therefore questionable. 

The data quality does not permit spectroscopy of individual 
flares. To characterize the spectral behavior of EXO 2030+375 
during the flaring episode and compare it with the rest of the 
outburst, we explored the luminosity-dependence of the source 
spectrum during the rising (flaring) part and the decay of the 
outburst. We grouped the individual INTEGRAL pointings ac- 
cording to the measured flux in the 20-80 keV range. For each 
group we extracted and analyzed the X-ray spectrum using the 
spectral model described above. The statistics did not allow us 
to explore the dependence of each individual spectral parameter 
on flux. We therefore fixed «h and Efon to their averaged val- 
ues (11 X 10^^ cm"' and 30keV, respectively) and explored the 
photon index F as a function of flux. 

The resulting dependence is shown in Fig.|6]for the rise (left) 
and decay (right) of the outburst. Higher values of F correspond 
to a softer spectrum (see the model description above). The hori- 
zontal error bars represent the width of the flux bins. The vertical 
error bars indicate Icr-uncertainties. While the spectrum appar- 
ently tends to get harder with increasing flux during the rising 
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part of the outburst, it remains roughly constant during the de- 
cay. To quantify this behavior, we performed linear fits to the 
data points and calculated the slopes with the corresponding un- 
certainties. The slope is (-9.0 + 3.7) x lO"-' (cts/s)"' in the rising 
part and (2.8 + 4.2) x 10"-' (cts/s)"' in the decay of the outburst 
(uncertainties at Icr confidence level). 

5. Discussion 

A direct comparison of the flares presented here with those 
observed with EXOSAT in 1985 revealed significant similar- 
ity. The peak fluxes and the relative amplitude of the flares are 
roughly the same in both episodes. Also the average shape - fast 
rise / slow decay - is similar in the two cases. The mean period 
is, however, different: ~4 hr for the EXOSAT flares and ~7 hr 
for the INTEGRAL ones. The location of the EXOSAT flares 
with respect to the nearest "normal" outburst is difficult to re- 
construct because no monitoring of the source flux (apart from 
the EXOSAT data themselves) was performed at the time. Since 
the o rbital phase of "nor mal" outbursts varies significantly with 
time dWilson et al.ll2002l) . the extrapolation of the orbital phase 
ephemeris back to the EXOSAT observations would not resolve 
the problem. 

X-ray flares in accreting pulsars are usually attributed to 
one of the following mechanisms: (1) instabilities of the ac- 
cretion flow around/within the magneto spheric boundary (e.g. 
iMoon et al.lF2003bl: iPostnov et alj|2008l) . (2) highly inhomoge- 
neous stellar wind of the donor star (e.g. iTaam et akl Il988t 
[Walter & Zurita Heras 200' ^, and (3) nuclear burning at the 
neutron star surface (e.g. L evine et al.l 120001 : iBrown & BildsterJ 
11998 ). The nuclear burning scenario is, however, difficult to 
reconcile with a relatively high accretion rate before and af- 
ter a flare, which would suppress the thermonuclear instability 
(Bildsten & Brown 1997). In the case of EXO 2030+375, inho- 
mogeneities of the companion's stellar wind are also unlikely to 
be a direct cause of the flares for the following reasons. First, the 
viscous time of the accretion disk that is be lieved to be presen t 
in EXO 2030+375 during normal outbursts (Wilson et al. 200% 
and t o even survive during quiescence (Havasaki & Okazaki) 
l2006h . is at least several days, which would smooth out any 
variations in the mass accretion rate M caused by inhomogene- 
ity of the wind shorter than this time. Second, nonuniform stel- 
lar wind cannot explain the observed quasi-periodic appearance 
of the flares. On the other hand, various kinds of magneto- 
hydrodynamic instabilities at the inner edge of the accretion disk 
may easily lead to oscillations in the mass flo w towards the po- 
lar caps of the neutr on star (Apparao 1991; Pos tnov et alJl2008t 
iD'A ngelo & Spruit 2010). leading to the observed flaring activ- 
ity. For example, D 'Angelo & Spruit (.201 0) have illustrated that 
when the magnetospheric radius (where magnetic field of the 
neutron star truncates the accretion disk) is larger but close to the 
corotation radius (where the Keplerian frequency is equal to 
the accretor's spin frequency), matter in the inner regions of the 
disk will pile up leading to an increase in the local gas pressure 
and, therefore, a decrease in r,„. When r,„ crosses r^, the accu- 
mulated reservoir of gas is accreted by the neutron star and the 
cycle repeats. 

To assess the applicability of this scenario to the flares ob- 
served in EXO 2030+375, we estimated the expected character- 
istic time scale of M-oscillations. Without going into physical 
details of the disk-field coupling at r,„ one would generally ex- 
pect that the oscillations in M take place on the time scale close 
to the local viscous ti me at ~ rc. Following t he standard a- 
viscosity prescription (IShakura & Sunvaevlll973h . this time can 



be estimated as ~ r^/v(rc) ~ l/[Qa(i//7?)^], where v(rc) is 
the viscosity at r^, Q is the spin frequency of the neutron star, 
and (H/R) the semithickness of the accretion disk. Using "stan- 
dard" values of a = 0.1, H/R = 0.05, and the known pulse 
period P ^ 40 s, one gets Tc ~ 7hr, i.e. of several hours, as was 
observed. 

The averaged profile of the flares is characterized by a steep 
rise and a slower decay (Fig.|4|, which is very s imilar to the 
flares observed with EXOSAT in 1985 (Fig. 2 of lParmar et"al] 
|J989). According to the authors, such a shape suggests a "drain- 
ing reservoir" that is in line with the picture described above 
(matter piling up on the inner edge of the disk). 

The observed difference in the spectral behavior between the 
flaring part and the rest of the outburst (Sect.|4]i suggests dif- 
ferent configurations of the region where matter couples to the 
field lines. Such a difference is indeed expected if the flares are 
caused by the oscillating inner edge of the accretion disk. In this 
case, matter from the oscillating inner disk rim would couple to 
different dipole field lines of the neutron star (and follow them) 
compared to the decay part of the outburst where the configura- 
tion of the inner disk rim is presumably stable. 

The difference in the mean period of flares in the 
INTEGRAL and EXOSAT observations can also be understood 
in the described picture. The period must depend on the mass 
transfer rate through the accretion disk, i.e. time needed to refill 
the reservoir. This rate could be different between the EXOSAT 
and INTEGRAL observations due to, e. g., changes in the state 
of the Be-disk. In Fig. |2] one might also notice some shortening 
of the the flare separation time as the flux increases (although 
this behavior is difficult to quantify with the available statistics). 
Such behavior, if real, might reflect shortening of the reservoir 
refill time as the mass transfer rate increases towards the maxi- 
mum of the outburst. 

Thus, we argue that the observational appearance of the 
flares in EXO 2030+375 suggests that the instability of the in- 
ner disk edge (pile-up/draining of matter) is the most probable 
cause of the flares. 

It is important to note that the rarity of the detected flaring 
episodes (even considering the relatively sparse observational 
coverage) means that the range of physical conditions needed 
to initiate flares could be very narrow, which would lead to the 
serendipitous character of the phenomenon. 
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